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We study the quantum Jarzynski relation for driven quantum models embedded in various envi- 
ronments. We do so by generalizing a proof presented by Mukamel [Phys. Rev. Lett. 90, 170604 
(2003)] for closed quantum systems. In this way, we are able to prove that the Jarzynski relation 
also holds for a bipartite system with microcanonical coupling. Furthermore, we show that, under 
the assumption that the interaction energy remains constant during the whole process, the relation 
is valid even for canonical coupling. The same follows for open quantum systems at high initial 
temperatures up to third order of the inverse temperature. Our analytical study is complemented 
by a numerical investigation of a special model system. 



I. INTRODUCTION 

If a system is driven far away from thermal equilib- 
rium by an external force, its behavior can no longer 
be described by linear response theory or other near- 
equilibrium approximations. Fluctuations may dominate 
the evolution of the system. Recently, some astonish- 
ingly general theorems, which make exact statements, 
have been found [H, 0, 11 i, 11 • For classically described 
systems there exist several proofs for various kinds of 
Hamilton dynamics [3, Q . 

One faces a different situation for systems the dy- 
namics of which are described by quantum mechanics. 
Mostly for very restricted situations or only for very spe- 
cial model systems have those theorems been shown to 
hold |i,0,[i,ll,[l3,llllll so far. However, there exists a 
general proof of the Jarzynski relation, a non-equilibrium 
work theorem, for closed quantum systems [13i] . Gener- 
alizations to open quantum systems have only been dis- 
cussed for restricted model systems or various constraints 
on the microscopic dynamics fill . [l3| . 

One major problem is that the proofs of the so-called 
group of non-equilibrium fluctuation theorems as well as 
the work fluctuation theorems relate properties of tra- 
jectories from forward and backward processes with each 
other. Trying to transfer these proofs to quantum sys- 
tems brings up the problem of trajectories in the quan- 
tum case. 

Another, possibly even more severe problem is the 
question of what exactly has to be understood by fluctu- 
ations in a quantum system. In quantum mechanics the 
state of a system is fully described by its density opera- 
tor. Without any measurement it is ambiguous to speak 
about fluctuations, because the probability of a measure- 
ment outcome of a certain quantity might vary, but not 
the untouched quantity itself. Fluctuations of one-time 
quantities are not ambiguous as far as their measure- 
ment is concerned. If, however, the relevant observable 
of two-time quantities, like work, e. g., does not commute 
with itself at different times then even the measurement 
of fluctuation becomes ambiguous. Every measurement 
changes the state of the system thereby influencing its 
evolution. In the case of work fluctuation theorems the 



measurement has naturally to be included since one needs 
to measure the amount of work performed on the system. 
This measurement will then yield different values for each 
repetition and is thus said to fluctuate. So, in this sense 
we have no problem speaking of fluctuations. 

This is why we try the following ansatz: We concen- 
trate on a work fluctuation theorem. In order to circum- 
vent the problem of the deflnition of a quantum trajec- 
tory, we use the quite general proof for closed quantum 
systems [l^, and try to generalize it to open quantum 
systems without having to define any quantum trajec- 
tory (in the sense of stochastic unravelling, i. e. continu- 
ous measurement). This will be done by splitting up the 
exponential function of the partition sum into its system 
and environmental part, respectively. There is no need 
to specify the underlying dynamics any further. It is not 
necessary to measure the system during the process since 
we rely on measurements at the beginning and the end 
of the process only. The idea of a two-time measurement 
approach has first been developed by Kurchan |6| . 

First, we start by a brief review of the Jarzynski re- 
lation in section ini Then, in section [HT] we will discuss 
conditions under which we will prove that the Jarzyn- 
ski relation is valid. This analytical study will then be 
supplemented by a discussion of our numerical results in 
section HVl 



II. JARZYNSKI RELATION 

If we perform a process on a system initially in a canon- 
ical state, the system will, in general, be driven out of 
equilibrium. Therefore, the average work W we have to 
bring up in order to perform the process will exceed the 
free energy difference AF between initial and final state: 



W > AF. 



(1) 



In 1997, C. Jarzynski came up with a remarkable relation 
that connects the non-equilibrium variable W with an 
equilibrium property of the system, AF [l[: 
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FIG. 1: Microcanonical coupling: between a small subsystem 
S and its environment. Total system is a closed compound 

system, initially in the canonical state ^(0) = 



Z(0) 



This relation holds no matter how far the system is driven 
out of equilibrium. It allows to access the free energy dif- 
ference between two states via non-equilibrium measure- 
ments. This is very useful for experirnents which cannot 
be carried out in quasi-equilibrium [l5| : Note that the 
right-hand side does not depend on how we get from the 
initial to the final state. 

The Jarzynski relation has been generalized to closed 
quantum systems by S. Mukamel [l3| under the condi- 
tion, that the system is initially in a canonical state. 



III. ANALYTICAL RESULTS 

We consider a bipartite system, split into the so-called 
system and environment, respectively. An arbitrary pro- 
cess is performed on the system leading to an explicitly 
time-dependent Hamilton operator. The compound sys- 
tem is considered closed and the complete Hamilton op- 
erator reads 



H{t) = H^{t) ® 



(3) 



where H^{t) denotes the Hamiltonian of the system, H'~^ 
the operator for the environment and if™* the interaction 
Hamiltonian with the coupling constant a, determining 
the strength of the coupling. In the following, we inves- 
tigate different coupling scenarios. 



A. Microcanonical Coupling 

We consider a microcanonical coupling [l^ (cf. Fig. [T|) 
first. This means that neither particle nor energy ex- 
change between the system and the environment is al- 
lowed. Note that because of entanglement between the 
system and the environment the work performed on the 
system will, nevertheless, be infiuenced by the environ- 
ment. This is due to the fact that the off-diagonal ele- 
ments of the reduced density operator of the system will 
be damped away. Since these elements are crucial for 
transitions within the system, the work should, indeed, 
be infiuenced by the microcanonical coupling. 



Since no energy transfer between system and environ- 
ment is allowed, we require that at any time t the Hamil- 
ton operator of the environment and the system, respec- 
tively, have to be separate constants of motion: 



Using [H^{t)^ H'-^] = we readily arrive at: 



0, = 0. 



(4) 



(5) 



The quantum Jarzynski relation for closed systems states 
that [cf. (HI)] 



-f3W 



m 
m 



(6) 



holds for arbitrary processes. Our aim is to prove that 
this relation also holds for this type of open quantum 
systems, i.e., for microcanonical coupling. To this end, 
we investigate the relation of the partition sums: 



Z(t) = Tr 



= Tr 



Here we have introduced the abbreviation = crT?'"*. 
Since the Hamilton operators commute with each other 
according to ([5]), we can split the exponential function, 
arriving at 



Z(i) = Tr 



(8) 



Using the definition of a canonical density operator 



z 

Z{t) = Tr 



we rewrite 



c ^-PH' 



(9) 



For the time-dependency of the trace we find based on 
the Liouville von-Neumann equation 



Z{t) 



dt [zS{t)ZC 
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Making use of the cyclic property of the trace and the 
commutator relations ([5]), we end up with 



At 



0^ 

ccan 



= ovt. 



(11) 



Thus, we particularly have 



Tr 



can^ 



Tr 



^fan(O) ® ^c„„e 



5^ . 
-can 
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(12) 



From this we conclude using Eq. ^ in Eq. 
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(13) 



The work on the compound system equals that on the 
open one since no energy exchange between the system 
and the environment is allowed and therefore the substi- 
tution in Eq. on the left-hand side is justified. Thus 
the Jarzynski equation holds for microcanonical coupling. 
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B. Canonical Coupling - Constant Interaction 
Energy 

Next, we turn to systems which are allowed to ex- 
change energy with their surroundings. The Hamilton 
operator is again of the form ([3]) . We assume the inter- 
action energy to be a constant of motion: 



H{t),H^ =0 



h^,h' =0. 



(14) 



This assumption is motivated by the fact that, for a suf- 
ficiently clear separation between system and environ- 
ment, the interaction energy should be very small com- 
pared to the energy of the system and the environment, 
respectively. Therefore, its change should also be consid- 
ered to be negligible. The proof then runs in complete 
analogue to that of the microcanonical coupling scenario. 
We start using the Jarzynski relation ^ for closed sys- 
tems. The partition sum can be written as in Eq. 
using the commutator relation ([1] 



(15) 



Z{t) = Tr[e~'5(^+^)1 # Tr[e-^-^e-'^^l , 



with A = H'^{t) + 1'^ ® and 13 = . Using 

[A, 13] — and re-substituting A and 13 we easily arrive 
at 



Z{t) - Tr ^l„(i) ® gg^e 



Z^{t)Z^ 



(16) 



as in the case of microcanonical coupling. Again, we 
investigate the time dependency of the trace using the 
Liouville von-Neumann equation and the cyclic property 
of the trace. 



Z{t) 



ZS{t)ZC 



= -iTr 



trcan 



{t) ® q\ 
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d 
dt 

(17) 

Using the commutator relation (fT4l) we conclude that 
[e-l^^' ,H(t')] = V i'. Thus, the relation ^ also 
holds for this scenario and therefore we conclude using 
Eq. (flB)) and relation pT|) in Eq. ([2]) that we still have 



-pw 



zs{Qy 



(18) 



This proves that the Jarzynski relation also holds for 
canonical coupling which conserves the interaction en- 
ergy. 

Note that the work APF performed on the compound 
closed system equals the work AWgub done on the sub- 
system S, i. e. AW = AVFsub [i3|- This is, in general, dif- 
ferent from the energy change of the subsystem during 
this process since the subsystem also exchanges energy 
with its environment 



sub 



(19) 



whereas the energy change of the total, closed system 
still equals the work performed on it: 



A&, 



AVF 



(20) 



So, the relation (jlSp is applicable if one would, in an 
experiment, measure ISW = AWgub, which is the work 
one has to supply in order to perform a chosen process. 
So, Eq. can be rewritten as 



z^io)- 



(21) 



C. Canonical Coupling - High temperature limit 

Finally, we consider systems with a relatively high tem- 
perature T, i. e. we have a small inverse temperature P 
at the beginning. We can then expand the exponent 
of the partition sum according to the Baker-Campbell- 
Hausdorff formula and related identities. In the follow- 
ing, we will use the following abbreviations: 



A = H'^it) (g) i"^ -f ® ijC" and B = aH"'\ 
We thus have 

which can be rewritten as 



^-PA^-PB ^ ^-P{A+B)^^[A,B]_ 



On the other hand, 

e-PB^^PA ^ ^-piB+A)^^^[B,A] 0(^3) 

We combine Eq. ([M]) with Eq. (^5]) and using 

-p{B+A) ^4lB,A] 



+ o{p') 



(22) 



(23) 



(24) 



(25) 



(26) 



we arrive at 



^-PA^-PB^-PB^-PA ^ ^^p{A+B)^lL-[A,B]^^[B,A]^-p{B+A) 



^ ^-2P(A+B) 



(27) 



Thus, we find: 



Z{t) = Tr 



,-PiA+B) 



'rB-!^B-!^A 



Tr[e-''^e-'3^] , (28) 



or, re-substituting A and B: 



Z{t) = Tr 



(29) 



We now can proceed as above: 

Z{t) = Tr[^^(t) ® fe^f^^'] Z^[t)Z^. (30) 
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FIG. 2: Model system for our numerical investigations. On 
the left-hand side the quantum well with a movable wall is 
depicted. The well is coupled to an environment with a single 
energy level of high degeneracy. 



We again investigate the time-dependency of the trace 
using the Liouville von-Neuman equation: 



d_ 
di 



Z{t) 



ZS{t)ZC 
= -iTr|[iJ^(i') 

= -iTr{[H^(t') 



'} 



(31) 



where we have used the cyclic property of the trace in 
the first line in order to get rid of the interaction term of 

the Hamiltonian. If we have H^{t'),H^{t) = V i', as 

for a time-dependent Zeeman splitting of a spin, e. g., we 
immediately arrive at 



d_ 

dt 



Z{t) 



ZS{t)ZC 



= 



and therefore 

Z{t) 



Z{0) ZS{Q) 



(32) 



(33) 



So, for high initial temperatures of the open quantum 
system, the Jarzynski relation holds if the Hamiltonian 
of the system commutes at different times. This can be 
shown to hold also in third order of f3, see appendix [X] 

From the derivation of this result one might expect 
that the Jarzynski relation would not hold in this form 
for arbitrary (3, but that it could fail and that the devi- 
ation of the average e~^^ from the relation of the par- 
tition sums should become larger the smaller the inverse 
temperature /3. This expectation will be tested numeri- 
cally. 



IV. NUMERICAL RESULTS 

As a simple model system we choose the standard po- 
tential well with width L containing one particle (particle 
in a box). Its spectrum is truncated at level n (cut-off). 



The well is coupled microcanonically to an environment 
(cf . Fig. ^ . To ensure that no energy is exchanged be- 
tween system and environment, we model the surround- 
ings as a single, highly degenerate, energy level 16]. 

One wall of this well (the left one, say) is now supposed 
to be movable, allowing for dilation and compression of 
the well [well width L{t)]. This time-dependence will 
model the working process. The particle is prepared ini- 
tially in a canonical state. Then the time-evolution of the 
system is calculated numerically under pure Schrodinger 
dynamics. From this we can infer the work distribution 
for any control function L{t). The discrete values of W 
result from the discrete spectral energy differences. The 
non-monotonous behavior of P{W) is a combined effect 
of the spectrum and the initial occupation probabilities. 
We use this P{W) in order to compute the average e~'^^. 
This result is then compared with the relation of the par- 
(t) 

tition sums -g^sjoj' "^l^ich can easily be computed indepen- 



dently for this model. We define e 



-f3W 



being a test factor. If the Jarzynski relation is correct, 
the test factor should always equal 1 within numerical 
accuracy. 

Since we want to consider the work done on the com- 
pound, closed system, we define the work as IF = AE, 
since the work simply equals the energy change of the 
closed system. Since we will, in general, find the system 
in different states for each measurement, the work is said 
to fluctuate. However, the definition of work and the 
corresponding fiuctuations is still a matter under dispute 



A. Closed Systems 

For our investigation of the closed system we choose 
the interaction strength with the environment to be zero. 

Our numerics show that, within numerical errors, for 
any process realization and any process velocity chosen, 
we, indeed, have 7j = 1 which means that the Jarzynski 
relation is found to be valid. This result holds for differ- 
ent wall displacement velocities, different initial temper- 
atures or different well widths. See Fig. [3] for a typical 
result for the work distribution obtained. 



B. Microcanonical Coupling 

Next, we turn to systems with a =/= 0. In order to 
check whether the Jarzynski relation holds, we need to 
consider a process for which the final Hamiltonian differs 
from the initial one; otherwise one trivially would have 



1. For now, wc focus on expansions of 



z{t) _ z'-it) 

Z{0) ~ 

the quantum well. 

We find that the work distribution is slightly affected 
by the microcanonical environment (cf. Fig. [¥]). For bet- 
ter comparison the work distribution for the closed sys- 
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FIG. 3: Distribution of work P(W) for a closed system. A 
circular process with -ff(tflnai) ~ H{0) and a maximal width 
of I/max ~ 21/(0) was chosen. The wall has been moved with 
constant velocity itv. The velocity of the process as well as 
the unit energy is in arbitrary units. Probabilities smaller 
than 1 x 10~^ are suppressed; cutoff at n = 5 is justified by a 
low initial temperature. The average work W is larger than 
zero. There is a finite chance to gain work while performing 
this process, as can be seen from the points to the left of the 
dashed line W = 0. 



FIG. 4: Distribution of work P(W) for a closed and a micro- 
canonical system, respectively. An expansion with i(tfinai) = 
21/ (0) was chosen. The wall has been moved with constant 
velocity v. The degeneracy was A'^'^ — 150 and the coupling 
strength a = 0.015 for the microcanonical coupling. For the 
closed system, a was set to zero. The velocity of the process 
as well as the unit energy is in arbitrary units. Probabilities 
smaller than 1 x 10"''' are suppressed; cutofi' at n = 5. 



C. Canonical Coupling 



tern is depicted, too. We get 7j = 1 within our numerical 
accuracy, indicating that the Jarzynski relation is valid. 

This also holds when changing the coupling strength, 
enlarging the environment by increasing the degeneracy, 
or changing the initial temperature. Note that the dis- 
tribution of work for the closed system is, indeed, though 
only slightly, different from the one obtained for the same 
boundary conditions under microcanonical coupling. 

Due to the interaction between system and environ- 
ment the energy levels of the system are " smeared out" . 
Despite the fact that the interaction is weak, the sys- 
tem's energy levels are each splitted into many energy 
levels with slightly different energy eigenvalues accord- 
ing to the number of energy levels in the environment 
and their respective degeneracies. Some of these energy 
levels are resolved within numerical accuracy leading to 
slightly distinct values of work W. This might lead to the 
impression that we have multiple values of work probabil- 
ities P{W) for a given work W in Fig.|4]to Fig. [7]for the 
open system. This is not the case since for exactly - in 
the numerical sense - equal work values, the probabilities 
add up to a single probability P{W). 



We finally turn to a canonical bath, investigating 
whether the Jarzynski relation still holds in the high 
temperature limit, as predicted. To this end, we cou- 
ple a canonical bath to the quantum well. It is modeled 
in such a way, that it allows every possible transition 
in our quantum well when in resonance [l6j . Its degen- 
eracy grows with increasing energy. For now, we focus 
on expansions of the quantum well, which is effectively 
cooled down by these dilations. This will cause energy 
to flow from the bath into the system trying to restore 
the initial temperature and we expect that transitions 
to higher energy levels will be facilitated. First, we have 
chosen a moderate initial temperature, such that only the 
first two energy levels are significantly occupied, resulting 
typically in a work distribution as depicted in Fig. [SJ 

We find that Eq. does not hold, the Jarzynski 
relation of this form is not strictly valid. We have for our 
test factor = 0.999, which quantifies the still small 
deviation from the Jarzynski relation. From sect. IIII CI 
we suppose that the lower the temperature the bigger this 
deviation. Indeed, if we choose the initial temperature 
to be very low, such that almost only the ground state 
is occupied, we have 7j = 0.98. Moreover, we expect the 
work distribution to be very narrow since transitions from 
the ground state into higher ones are relatively unlikely 
(cf. Fig.©. 

Next, we choose a very high initial temperature, such 
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FIG. 5: Distribution of work P(W) for a canonical system. 
An expansion with I/(tenai) = 2L(0) and a = 0.015 was cho- 
sen. The wall has been moved with constant velocity v. The 
velocity of the process as well as the unit energy is in arbi- 
trary units. Probabilities smaller than 1x10"^ are suppressed; 
cutoff at n = 5. 



FIG. 7: Same as Fig [5] but with high initial temperature 



V. CONCLUSION 



0.1 
0.01 

E 

0.001 

1x10"'' 
1x10-5 



■ cold * 




X 




X 




X 




X 




X 




X 




* XX 




* X 

X S X 
1 Y 1 1 


1 



Work W 



Starting from the quantum version of the Jarzynski re- 
lation for closed systems, we have shown that this relation 
can be generalized to open quantum systems of different 
coupling types. The key idea has been to split the expo- 
nent of the partition sum into parts for the system and 
environment, respectively. This is, in general, not possi- 
ble, since the system and environmental Hamiltonian do 
not commute with the interaction part. This fact would 
result in additional terms containing the commutator of 
these Hamiltonians. 

Nevertheless, we have been able to split up the par- 
tition function under the selected conditions, since then 
the correction terms cancel. 

We have stressed that the problem of fluctuation is 
intimately related to some measurement scheme. Further 
investigations should be guided by the observation that 
in quantum mechanics it seems more natural to consider 
the fluctuation of a (subsystem-)energy rather than of 
work proper . 



FIG. 6; Same as Fig. O but with low initial temperature 

APPENDIX A: THIRD ORDER EXPANSION IN 
TERMS OF /3 < 1 



that almost every energy level has the same probability 
to be occupied, in order to check whether the Jarzynski 
relation is valid as predicted. The work distribution is, as 
expected, rather broad for these temperatures, cf. Fig.[7l 
Nevertheless, we find that 7j = 1 within numerical accu- 
racy. This means that Eq. is fulfilled for high initial 
temperatures, confirming our high-temperature limit. 



We want to extend the result of sect. IIII Cl to the third 
power of We proceed analogously as above, trying 
to split up the exponential function. For 0(/3*) = we 
arrive at 

g-/3Ag-/3B ^ ^-f}(A+B) + SllA,B]^-!^{[A^lA,B]] + [B,lA.B]]) 

(Al) 
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and 



and analogously 



-PB -PA ^ p!^{[A\A,B\] + [B\A.,B\]) -p(B+A)-Sl[A,B] 



Combining Eq. (jAip with Eq. \A2\ we have 



(A2) 



g-/3C_g-/3C+ 



e^^^e-2^(^+^)+0(/34), (A8) 



^-I3A-13B-I3B-I3A _ I3{A+B) + S1[A,B]^-I3(B+A)-!^[A,B] 



ED ED 



g-/3C+g-/3C_ 



(A3) 



with C+ = -I3{A + B) + ^[A,B] and C_ ee + 

J3) — ^[A, S]. We investigate the resulting product in- 
troducing the abbreviation [C+,C_] = D±: 

^-P(C++C-) + ^D±-^{[C+,D±] + [C-,D±\) J^0{I3^). (A4) 

Since C+ = C_ + (3'^[A, B] we have 

D±^-P^[[A, [A,B]] + [BM,B]]] 



(A5) 
(A6) 



This leaves us with 



^-PC+^-pc. 



-P{C++C^) + I2^D± 



-2piA+B) ii^D± 



e V^± + 0(/3*) (A7) 



with Dzp = [C-, C+] = —D±. Combining these two rela- 
tions we have 



g-/3C+g-/3C_g-/3C_g-/3C+ 



Therefore we have 



-4P{A+B) 



(A9) 



Trfe-'^-4g-;3SN 



m / -l^A ~S.B -S.R -S.A 

= Tr e 2^e ^ e ^ e = 



Trfe~'^+e~'" 
Tr^eT^'^+e^'^ 

Tr(e-'^(^+^)) = 



e * e 4 



(AlO) 



The rest of the proof for the Jarzynski relation to hold 
runs analogously to sect. IIII Cl 
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